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We are interested in the existence of depolarization waves in the human brain. 
These waves propagate in the grey matter and are absorbed in the white matter. 
We consider a two-dimensional model ut = Au + f(u)l\ y \<R — aul^^R, with / a 
bistable nonlinearity taking effect only on the domain M. x [-R, R], which represents 
Q^ ■ the grey matter layer. We study the existence, the stability and the energy of non- 

-^ trivial asymptotic profiles of possible travelling fronts. For this purpose, we present 

dynamical systems technics and graphic criteria based on Sturm-Liouville theory 
and apply them to the above equation. This yields three different behaviours of the 
solution u after stimulation, depending of the thickness R of the grey matter. This 
may partly explain the difficulties to observe depolarization waves in the human 
brain and the failure of several therapeutic trials. 

Keywords: spreading depression, reaction-diffusion equation, travelling fronts, 
pg ' Sturm-Liouville theory. 

r- AMS classification codes (2000): 34C10, 35B35, 35K57, 92C20. 
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1 Introduction 

The propagation of depolarization waves, also called spreading depressions, may appear in 
a brain during strokes, migraines with aura or epilepsy. In rodent brains, the propagation 
of depolarization waves has been observed for more than fifty years [15]. During stroke 
in rodent brain, they cause important damages and are therefore a therapeutic target. 
Pharmacological agents blocking the appearance of those waves have been studied and 
reduce strongly after-effects of stroke in the rodent brain [TS1 [TH] . However propagation of 
depolarization waves during stroke in the human brain is still uncertain [TJ [161 H3 E31 EI] 
and the pharmalogical agents used in rodent have seemed to have no effect on human 
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stroke. Mathematical models of these depolarization waves may help to understand these 
points. 

In [2], the first author has build a mathematical model of these waves and has studied 
numerically how the morphology of the brain may influence their propagation. Simplifi- 
cations of the mathematical model lead to the study of propagation phenomena for the 
following biological equation 

— -Au = f(u)l n -aul RN \ n , teR,XeR N , (1.1) 

where f(u) = A u{u — a)(l — u) is the usual bistable nonlinearity with a g]0, l/2[ and 
A > 0, and where a is a positive number. We denote by In the caracteristic function of 
the domain Q that is 1q(X) — 1 if X G £1 and 1q(X) = elsewhere. The function u 
represents the depolarization of the brain that is that if u(X) = the brain is normally 
polarized at the point X whereas if u(X) = 1 the brain is totally depolarized. The domain 
Q represents the grey matter of a brain, where the react ion- diffusion process that triggers 
depolarization waves takes place, whereas M N \ Q contains the white matter, where the 
waves are absorbed. The problem is to understand the influence of the geometry of Q 
on the propagation of waves. In particular, the layer of grey matter in the human brain 
is thiner than in other species and admits more circumvolutions. The influence of the 
circumvolutions of the grey matter is partially studied in [I]. 

We will focus here on the influence of the thickness of the layer of grey matter. We 
may assume that Q is a straight cylinder of radius 2R. In [2], numerical studies have 
shown that small values of R may prevent the spreading of the ionic waves. A partial 
theoretical study in any dimension N has been led in [3] , where the first author has proved 
the non-existence of the depolarization waves if the thickness R is small enough and their 
existence if R is large enough. 

The results of [5] are not completely satisfactory. First, they only deal with R small 
or large enough. Secondly, no complete description of the possible asymptotic profiles of 
the travelling waves, their stability, or their energy, has been pursued. As a consequence, 
even if the existence of waves is proved, nothing is known about their asymptotic profiles 
and their stability. The purpose of the present paper is to complete this study for any R, 
in dimension N = 2. More precisely, we consider the equation 

— = Au + f(u)l ]y] < R -aut ]y]>R , teR,(x,y)eR 2 . (1.2) 

We are looking to solutions travelling in the x direction at speed c, that are solutions 
u of fll.2p which can be written u(x,y,t) = v(x — ct,y). Travelling fronts are solu- 
tions u(x,y,t) = v(x — ct,y) such that there are two asymptotic profiles VI ^ V + with 
lim^-too ||i>(£, .) — V±||hi(h.) = 0. Using standard elliptic estimates, a profile V is a solution 
in H 1 (M) of the equation 

V" + f(V)l lyl < R - aVl lyl>R = 0, y G R . (1.3) 

Notice that a profile is trivially associated to an equilibrium point E(x,y) = V(y) that 
is a stationnary solution of (11.21) . The trivial equilibrium point E = corresponds to 
the normal state of the brain, whereas a non-trivial equilibrium point corresponds to a 



depolarized state, deleterious for the brain. We are more particularly interested in the 
existence of travelling fronts with positive speed c connecting the profile V+ = with a 
non-trivial profile V-. Such fronts correspond to the invasion of the equilibrium state 
by a deleterious state. 

The purpose of this article is to achieve a complete study of the existence of the 
profiles, their stability and their energy (see the following sections for the definitions of 
stability and energy). 

Main result 1.1. There exist two critical thicknesses Ri > Rq > such that: 

i) if < R < Rq, there is no non-trivial profile, i.e. that V = is the only solution of 

(Ojl inH x (R). 

ii) If Rq < R, there exist non-trivial profiles. One of them, denoted by Vm, is larger than 

every other one. The largest profile Vm is stable, and every other non-trivial profiles are 

unstable. 

Hi) The energy of the unstable profiles is always larger than the energy of the stable profiles 

and Vm- If Rq < R < R\, the energy o}Vm is larger than the energy ofO, whereas it is 

smaller if R> R\. 

In this paper, we first adapt several graphic criteria for stability of profiles, based 
on Sturm-Liouville theory. Then, we present dynamical systems technics for studying 
energy of profiles. These criteria and technics are general ones and are valid for any 
nonlinearity /. We obtain the above results by applying these ideas to the particular 
nonlinearity f(u) = X u{u — a)(l — u). The check on the graphic criteria sometimes relies 
on a numerical study of the phase plane. 



Once the profiles and their stability completely described by the Main Result 11.11 one 
is able to obtain the existence of travelling fronts. 

Consequences 1.2. 

i) If < R < Rq then there is no travelling front for Equation fll . 2j) . 

ii) If Rq < R < R\ there exists a globally stable travelling front u(x,y,t) = v{x — ct,y) 

with rini£_> +00 v(£, .) — and lim^^oo v(£, .) = Vm with a negative speed c < 0. 

ii) If Ri < R there exists a globally stable travelling front u(x,y,t) = v(x — ct,y) with 

lini£_> +00 v(£, -) = and lim^.oo v(£, .) = Vm with a positive speed c > 0. 



Assertion i) is a clear consequence of the Main Result 11.11 The proof that this Main 
Result also implies the other assertions is not the goal of this paper and will not be 
detailled. One way to prove it is to use a method recently introduced by Risler in [20] . 
The fundamental idea is to use the existence of an energy functional for (11.21) in any 
Galilean frame travelling at constant speed s in the x direction. Keeping this basic idea, 
the original proof has been revised in [TD] and [9]. The adaptation of these technics to 
Equation (\1.2\i has been detailled in [3] . The original motivation of Risler's technics was 
to develop a method of proof of existence and stability of fronts in equations where no 
comparison principle is available. As suggested by the work of the first author in [3], the 
method of Risler is also usefull to get the existence and stability of fronts in equations as 
(jl.2p for which even if there is a comparison principle, no particular positive subsolution 
is known. 



The results of this article yield a good insight into the different behaviours of the 
depolarization in the grey matter after it has been stimulated. If the stimulation takes 
place in a part of a brain where the grey matter is thin (R < Rq), the neurons of the 
grey matter quickly repolarize, i.e. the solution of ( II. 2p goes back to zero, uniformly and 
exponentially fast. If the grey matter is slightly thicker (R < R < Ri), the repolarization 
is slower and more progressive, the depolarized area disappearing by shrinking. Math- 
ematically speaking, there exists a stable non-trivial profile but with an energy larger 
than zero. Thus, the equilibrium state zero invades the excited state by travelling fronts, 
reducing the excited area. Notice that it may take very long time to go back to rest if the 
initial excited area is large. Finally, if the grey matter is thick (Ri < R), a depolarization 
wave propagates. Indeed, the excited state has a lower energy than zero. These three 
behaviours are illustrated in Figure [TJ The fact that the behaviour depends on the width 
of the grey matter may explain why the depolarization waves have been observed or not 
in the human brain depending on the experiments, see [H[T6l[T7] or[23] and the discussion 
of Section 





Figure 1: The three typical behaviours of the depolarization in the grey matter after stim- 
ulation. From left to right: R < Rq, Rq < R < R± and R\ < R. 



The paper is splitted as follows. In Section [2] we explicit a relation between the profiles 
and the equilibrium points of a parabolic equation in (—R,R). In Section [31 we obtain 
graphic criteria of existence and stability of the profiles, mainly based on Sturm-Liouville 
theory. The energy of the different profiles is studied in Section HI using technics coming 
from the infinite dimensional dynamical systems theory. Finally, in Section \5\ we discuss 
the relations between our mathematical results and the biological phenomena of spreading 
waves in the brains. 

2 Relations with a parabolic equation on (— R, R) 

In this preliminary section, we enhance an obvious relation between the profiles and the 
equilibrium points of a parabolic equation on the segment (—R,R). We also recall some 
basic properties of the dynamics of this parabolic equation. In fact, in this paper, we 
could perform all the arguments with the parabolic equation on the whole domain K. 
However, this would bring useless technicalities. 

2.1 Back to a bounded interval 

The profiles V are the solutions of Equation (11.31) . In other word, they are the stationnary 
solutions of the evolution equation 



du d 2 u 
dt dy 



. g{v,u) (2.1) 
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where g(y, u) = f(u)l\ y \<R— cral| y |>.R. The energy corresponding to this reaction-diffusion 
equation is the functional E : iJ 1 (M) — > R defined by 



E(u) 



du 
dy 



(y) 



dy 



G(y,u(y))dy 



where G(y, u) = j Q g(y, v)dv is a primitive of g. 



(2.2) 



We notice that any profile V satisfies V(y) = c + e^ v + c_e _v/ " y outside (—R,R). 
Since V belongs to H 2 (M.), one has c + = for y > R and c_ = for y < R. Since 
iJ 2 (R) C C 1 (R), V is continuous and solving (jl.3p is equivalent to solving the following 
problem 

' ^(2/) + /(y(2/)) = ye(-R,R) 

< V'(-R) = ^V(-R) (2.3) 

k V'(R) = -v^ y(i2) . 

Equation (12. 3p caracterises the equilibrium points of a parabolic equation on (-R, R) 
with Robin boundary conditions 



du 

~dt 
du 

dy 
du 

I dy 



(y,t) 



d 2 u 



dy 



2 {y,t) + f(u{y,t)) (y,t)e(-R,R)x 



(-R,t) 

(R,t) = 



a u(—R,t) 
a u(R,t) 



t> 
t > 



(2.4) 



The energy corresponding to Equation (12.41) is given by 



S(u) 



R 



du 

dy 



(y) 



dy 



R 



F(u(y))dy+^{\u(R)\ 2 +H-R)\ 2 ) 



(2.5) 



where F{u) = J " f(v)dv is a primitive of /. 

It is straightforward to verify that if V is a profile, then E(V) = £ (V^-j^ff^). As a 
conclusion, the problems of the existence and the energy of the profiles are equivalent to 
the study of the equilibrium points of the parabolic equation (12. 4p and their energy. We 
will often abusively denote V^r^) by simply V. However, notice that the dynamics of 
the equation (12.11) on the whole real line and the ones of (12. 4p on (-R, R) are different. 
Moreover, the spectrum of the linearization of (12. ip at a profile V is different from the 
spectrum of the linearization of (12 .4p at V\(-nm, even if we prove in Section I3T21 that the 
number of positive eigenvalues is the same. 



2.2 Basic facts about the parabolic equation on (-R, R) 

The theory of the parabolic equation (12 .4p with Robin boundary conditions is very well- 
known, see for example [14] . We recall here some basic dynamical properties, which will 
be used in this paper. 

First, (12.41) admits a comparison principle. 



Lemma 2.1. Let u(t) and v(t) be two solutions of (12. 4ft with respective initial data uq 
and vq belonging to H 1 (—R,R). Assume that uo(y) < vo(y) for all y G (—R,R). Then, 
for all (y,t) G (-R,R) x R + , u(y,t) < v(y,t). 

Secondly, the parabolic type of (12 .4p and the fact that f(u)u < for large u im- 
plies that any solution u(t) with initial data uq G H x {— R, R) is bounded in H k (—R,R) 
uniformly for t G [to, +oo) for any to > and for any k G N. This yields a compacity 
result. 

Lemma 2.2. Let (t n ) n£ N be a sequence of times such that liminft n > 0, let (w°)„eN C 
H l (—R, R) be a sequence of initial data and let u n (t) be the solutions of (12.41) with u n (0) = 

Then, there exist an extraction ip and a function u M G ^(—R, R) such that ||w^( n )(t^( n )) — 
Uoo\\ m{~R,H) — > when n goes to +oo. 

The energy £ defined by (12. 5p is a strict Lyapounov functional for the dynamics of 
(12. 4ft : for each u(t) G H l {—R,R) solution of (12. 4ft . the function £{u(t)) is decreasing in 
time, except if u(t) is constant, that is except if u is an equilibrium point, and thus a 
profile. Indeed, it is straightforward to get d t £(u(t)) = — J_ R \^(y,t)\ 2 dt. The existence 
of a strict Lyapounov function yields what is called gradient dynamics. In addition of the 
compacity property of Lemma [2.21 it is very classical to obtain Lasalle principle, see 



Lemma 2.3. Let u{t), t > 0, be a solution in iJ 1 (— /?, i?) of (12. 4p . Let u(u) be the 
uj— limit set of u, that is 



U(U) 



a\u) 



\ v G H l (-R, R) J 3(t n ) n6N , t n > +oo and \\u{t n ) - v\\hU-r,r) > \ 

Then, u(u) is a non-empty connected compact set which only consists in equilibrium 

points, i.e. in profiles. 

Similarly, if u{t) is defined for all t G IR and uniformly bounded, the a— limit set 

I v G ^(-R, R) J 3(t n ) neN , t n y -oo and \\u(t n ) - v\\ h ^-r,r) > \ 

is a non-empty connected compact set which only consists in profiles. Moreover, if V_ G 
a(u) and V + G u(u) then £(VJ) > £(V + ) and in particular V_ ^ V+. 

3 Existence and stability of profiles 

We are interested in the asymptotic profiles of the possible travelling fronts of Equation 
(II. 2p . We recall that a profile V G iJ 1 (R) is solution of (11.31) . We wonder if there exist 
non-trivial profiles (i.e. profiles V j£ 0) and if they are stable. A profile V is stable if the 
linear operator L v : H 2 (R) — y L 2 (R) defined by 

L v tp = </ + (f{V)l\ y \< R - al lyl>R ) up (3.1) 

has no spectrum in the half-plane {z G C, Re(z) > 0}. 

In this section, we give graphic criteria for existence and stability of profiles and we 
apply them to Equation ( II. 3p to obtain the following result. 
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Proposition 3.1. For any given a, X and a, there exists a positive number Rq such that: 

i) if R < Rq, then there does not exist any profile different from 0. 

ii) if R> Rq, then there exist non-trivial profiles. One of them, denoted by Vm is strictly 

larger than every others, smaller than 1 and is even. Moreover, if R > R Q then Vm and 

are stable profiles and are the only ones to be stable. 

Hi) if a > aX then if R> Rq, there exist exactly two non-trivial profiles. 



3.1 A graphic criterium for existence of profiles 

We recall here the simple graphic arguments already mentionned in [3]. As already no- 
ticed, looking for profiles V solutions of ( II. 3p is equivalent to looking for solutions V of the 
equation (12. 3p . In the phase plane, the flow $> yo : M 2 — > M? associated to the differential 
equation V" + f(V) = between —R and y is given by 



$ yo (V , Vi) = (V(y ), V'(y )) where V solves 



V"{y) + f(V{y)) = 0, ye(-R,y ) 
(V(-R),V'(-R)) = (V ,V 1 ) 



Let V and V C M 2 be the straight lines R.(l,y/a) and R.(l,-y/a). Notice that the 
conditions V'(-R) = ^a V(-R) and V'(R) = -y/a V(R) correspond to (V,V')(-R) G 
V and (V,V')(R) E V respectively. Hence, the graphic interpretation of ( 12. 3 \ is the 
following. 



Graphic criterium 3.2. The profiles, i.e. the solutions of f 1 2 . 3 1) . correspond to the 
trajectories of the flow $ y that start on T> when y = —R and finish on V when y = R. 
Therefore, they correspond to the intersections of the curve ^rD with the straight line T>' . 



The results of existence stated in Proposition 13.11 follow from Graphic Criterium 13.71 
The critical radius Rq is the first (positive) R such that QrD intersects T>' elsewhere than 
at (0,0). In Figures [2] and El we illustrate the phase plane, the flow Q y and the curve 
$>rD for different cases. Figure @] shows the graph of four profiles V(y) in [-R, R] and 
the corresponding intersections in the phase plane. 




Figure 2: The phase plane for a > aX and &rT> for three increasing values of R: two 
cases with no non-trivial profiles and the last with two non-trivial profiles. 



Figures [2] and [3] show that there are two cases. If T> is always outside the domain 
delimited by the homoclinic orbit to zero, then there cannot exist more than two non- 
trivial profiles. If T> intersects the homoclinic orbit to zero, then for R large enough 



Figure 3: The phase plane for a < aX and <&rD for three increasing values of R: no 
non-trivial profiles, two ones and four ones. 





Figure 4: The intersections of §rD and T>' in the phase plane and the corresponding 
profiles in the (y, V) plane 



there exist as much non-trivial profiles as wanted. The critical value of a is obtained 
by computing the unstable direction of in the phase plane. This direction is given by 

1 



the eigenvector (1, v aX) related to the eigenvalue v aX of the matrix 



aX 



which 



corresponds to the linearization of the flow at the equilibrium point (0,0). This shows 
that T> intersects the homoclinic orbit to zero if and only if a < aX. 

The fact that there exists a profile Vm larger than the other ones is classical, see for 
example [21J. It is a consequence of Lemmas 12. II and 12. 31 if u(t) is a solution of (12 .4p with 
u(0) larger than any profile, then the cu— limit set u(u(0)) consists in a profile Vm larger 
than any profile. The fact that Vm is even follows from the symmetry of the problem 
since Vm{— • ) is also a profile. Notice that this symmetry corresponds to the symmetry of 
the phase plane and to the fact that the trajectory y h> $ 2/ (Vm(— R), V m (—R)) intersect 
T>' for the first time at y = R. 



3.2 A graphic criterium for stability of the profiles 

Let V be a profile solution of (II .3p . We consider the linearization Ly of the flow along 
the profile V. We recall that Ly is given by (13. ip . The operator Ly is self-adjoint and is 
a relatively compact disturbance of A — a and therefore has the same essential spectrum, 



see [H]. This yields the following characterisation of the stability of a profile V, see [22] . 

Proposition 3.3. The spectrum of Ly consists in ] — oo, —a] and a finite number of 
isolated real eigenvalues pi > —a of finite multiplicity. As a consequence, the profile V is 
asymptotically stable if and only if Ly has no nonnegative eigenvalue. 

We are looking to nonnegative eigenvalues of Ly that is to A > such that there exists 
ip G H 2 (R) such that ip ^ and 

V"{y) + (f'(V(y))l mR - al lyl>R ) <p(y) = pp(y) • (3-2) 

By the same argument used for V in Section [2J <p> is explicit outside (— R, R) and it is 
equivalent to look to ip satisfying 

<p"(v) + f'(V(y)My) = Mv) v e {-R, R) 

p'(-R) = TjJTTmp(-R) (3.3) 

ip'(R) = -^JTT^p(R) 

In this section, we adapt the standard Sturm-Liouville theory to the eigenvalue problem 
(13. 3p . Keeping in mind the geometric interpretation of Sturm-Liouville arguments, we will 
obtain a graphic criterium to count the number of nonnegative eigenvalues of Ly. The 
idea of using such graphic arguments for studying the one-dimensional parabolic equation 
is not new, see [S]. In fact, one can understand the whole dynamics of a one-dimensional 
parabolic equation by similar arguments as shown in [B]. See for example [7] for a review 
on this subject. 

For any fi > we define # M by 



9^ = -sm 2 9, + (fi-f'(V(y)))cos 2 9^ , ye{-R,R) 
6^{—R) = arctan(- v //T+~a7) 



(3.4) 



and we set for fi > — a 



/i(/i) = 9^{R) + arctan( v / y u + a) . (3.5) 



Lemma 3.4. A nonnegative number ft is a nonnegative eigenvalue of Ly if and only if 
h(fi) belongs to 7rZ. 

Proof: Let [i > and ip G C 2 ([— R, R]), ip ^ 0, satisfying (I3.3p . We introduce two 
functions of class C 2 , p > and 9 G E, such that (<p,<p')(y) = p(y) (cos 9 (y), sin 9 (y)). 
Notice that these functions exist and are regular since the vector (</?, ip') cannot vanish 
because p is a non-trivial solution of a linear second-order differential equation. Up to 
the change of the sign of <p and to the addition of a constant in 7rZ to 9, p'(—R) = 
y/p + ap(—R) implies 9(—R) = arctan( A //7+~a7). Moreover, it follows from (I3.3P that, 
for y E(-R,R), 

9' = - sin 2 9 + (n- f{V{y))) cos 2 9 

p = p sinfl cos 9(l + n-f'(V(y))) (3.6) 

Hence, 9 = 9^ and p'{R) = —\/p + ctp{R) is equivalent to # M = — arctan(^//i + a) + kir 
with fceZ. 
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To show the other implication, it is sufficient to follow the previous arguments in the 
opposite way: we start with a function 9^, we construct p by (13. 6p and we check that 
ip = pcosdp is a solution of ( I3.3p . D 



Lemma 3.5. Let V be a profile solution of ( II. 3p . 

i)Ifp> ||/'(^)lloo then, for ally £[-RiR\, O^y) e {U,l) and h(p) e (0,vr). 

izj I//ii > /i2 i/ien, /or a// y G [— i?, i?] ; 9 tJ , 1 (y) > d^iv) an d h is an increasing function. 

Proof: We prove assertion i) by contradiction. Let yo = inf{y G [—R,R]/9 /J ,(y) $. 
(0,tt/2)}. Notice that y > -R since 9^-R) G (0,vr/2). If 9^y ) = 0, then 0^(y o ) = 
t 1 ~ f'(V(yo)) > which is absurd since 9^(y) > for y < y by definition of y . The 
contradiction is similar if 9^(y ) = n/2. 

The proof of ii) is very similar: let y = inf{y G [—R,R]/9 fJ , 1 (y) < 9^ 2 (y)}. We have 
0« (yo) = W (y ) . If W (y ) & n/2 + ttZ, then 9'^ (y ) - 9'^ (y ) = (pi - p 2 ) cos 2 W (y ) > 
and it contradicts the definition of yo. If O^Do) G 7r /2 + 7rZ, then 9' /xl (yo) = 0'^ 2 (yo), 
PpM = e '^(Vo) an d 0£(yo) = 2 + /n - f(V(y )) > 9 f ^(y ). In both cases, we obtain 
the desired contradiction. □ 



As a consequence of both preceding lemmas, we obtain the following Sturm-Liouville- 
type result. 

Proposition 3.6. Let V be a profile solution of ( II. 3p . Let Ly be the operator defined by 

( 13. ip and let h be the function defined by ( 13. 5p . 

i) k G N is such that h(0) G (—kit, (1 — k)ir] if and only if the operator Ly has k nonneg- 

ative eigenvalues. 

ii) if Ly has a nonnegative eigenvalue, there exists a positive eigenfunction associated to 

the largest eigenvalue of Ly. 

Hi) the nonnegative eigenvalues of Ly are simple. 



Proof: By Lemma 13. 4[ the nonnegative eigenvalues correspond to the values p for which 
h(p) G Z. Since, as shown in Lemma 13. 5| p H- h(p) is increasing and belongs to (0, tt) 
for large p, then Assertion i) follows from the continuity of h. Moreover, if it exists, 
the largest eigenvalue p is necessary such that h(p) = 0. Thus 9^(R) G (— 7r/2,7r/2). 
Since 9^(—R) G (—tt/2,tt/2) and since 9' (y) < as soon as 9^(y) = ±7r/2, we must 
have 9^(y) G (— it/2, it /2) for all y G [— R, R]. Following the arguments in the proof of 
Lemma 13.41 we construct an associated eigenfunction tp = p cos 9^ which is positive on 
[— R, R]. It is positive on R by the extension imposed by (13. 2p . Finally, the eigenspace 
of a nonnegative eigenvalue is one-dimensional since any eigenfunction must satisfy (13.31) . 
Since Ly is self-adjoint, any nonnegative eigenvalue is simple. □ 



The interest of the first assertion of Proposition 13.61 is to have a graphic interpretation. 
Indeed, let <p be the solution of (<p,<p')(—R) = (1, yfa) and ip" + f'(V)<p = 0. In other 
words, y H- (ip,ip')(y) is a trajectory of the linearization of the flow Dy§ y along the 
profile V such that (tp,tp')(—R) is in the tangent space of V at (V,V')(—R). For each 
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y G (-R, R), (<f, <p')(y) is therefore in the tangent space of <& y U at (V, V')(y). By the same 
arguments as in the proof of Lemma f3.4| 6o(y) defined by (13.41) is exactly the argument 
of the vector (ip,(p')(y). Therefore, one can observe the angle 6o(y) by looking at the 
argument of the vector tangent to <&yT> at (V, V')(y). Thus, we can graphically interpret 
Assertion i) of Proposition 13.61 in the following way. 



Graphic criterium 3.7. Let V be a profile, solution of (II. 3p and let Ly be the operator 
defined by (13. 1| . Let y G [-R, R] i — > Ay(y) G K 2 be a continuous function such that for 
all y G [— R, R], Ay(y) is a unit vector tangent to Q y T> at (V,V')(y). Then, the number 
of times Ay(y) crosses the direction ofV when y describes [—R,R] is exactly the number 
of nonnegative eigenvalues of Ly. The crossings have to be counted in a algebraic way: 
positively in the clockwise sense and negatively otherwise. 

The previous graphic criterium is illustrated in Figure [5j This criterium requires the 
knowledge of the curve <& y T> for each y G (—R,R). Therefore, it is difficult to apply as 
it stands. Fortunately, there is a way to count the number of times Ay(y) crosses the 
direction of T>' with the knowledge of the curve &rD only. Indeed, in the particular case 
of this article, the phase plane associated to / is such that the curve §rD is above T>' 
in {(u, v) G M 2 , u > 1} and below T>' in {(u,v) G M 2 , u < 0}. For this reason, using 
homotopy and topological arguments, the vectors A (y) and Ay M (y) corresponding to the 
extremal profiles and Vm must have a trivial number of crossings with the direction of 
V . Moreover, for another profile V, we can compute the number of crossings of Ay(y) 
with T>' as follows. 




Figure 5: an example of application of the graphic criterium \3.1\ One counts along 
a profile V the number of crossings between the tangent vector to the curve & y T> and 
the direction of V . This gives the number of nonnegative eigenvalues of Ly. Left, the 
largest profile has no nonnegative eigenvalue. Right, the other non-trivial profile has one 
nonnegative eigenvalue. 



Graphic criterium 3.8. Let V be a given profile. We follow the curve QrD from (or 
from Vm) to V . We count in a algebraic way the number of times the unit tangent vector 
to QrD crosses the direction ofT>' during this course. The resulting number is exactly the 
number of nonnegative eigenvalues corresponding to the profile V . 
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The equivalence of the graphic criteria 13.71 and 13.81 can be easily "seen on the figure" . 
However, the rigourous proof uses homotopies and the topology of the plane, similarly to 
Jordan theorem. This proof is not the subject of this paper and will be omitted. 

We can apply the above graphic criterium to the particular phase plane of this article. 
Because the flow of V" + f(V) is turning around the equilibrium (a, 0), every profile V, 
which is not one of the extremal profiles or Vm, is unstable, see Figure EJ 




Figure 6: examples of application of the graphic criterium \3.8l One follows the curve 
QrD from to a profile V and counts in a algebraic way the number of crossings between 
the tangent vector to the curve and the direction of T>' . One then obtains the number 
of nonnegative eigenvalues of Ly. At the left, the largest profile (2) has no nonnegative 
eigenvalue whereas the other non-trivial profile (1) has one nonnegative eigenvalue. At 
the right, the four non-trivial profiles have respectively one, two, one and zero nonnegative 
eigenvalues. 



4 Energy of the profiles 

We recall that the relevant energy for the profiles is the function E defined by (12.21) . We 
also recall that the energy E{V) of a profile V is equal to the energy ^V^-^)) of the 
restriction of V to (—R,R), see Section [2J Let Rq be the critical thickness introduced in 
Proposition 13.11 We know that for R > Rq there exists a unique non-trivial profile Vm 
which is stable and that Vm is the largest profile. This section is devoted to the following 
result. 

Proposition 4.1. 

i) IfV is an unstable profile, then E(V) is larger than E(0) = and E(Vm)- 

ii) There exists R% > -Ro such that on (Rq, R2) (this set being possibly empty, for example 

if a > aX) the energy E(Vm) of the largest profile is an increasing function of the radius 

R and on (R 2 ,+oo), E(Vm) is a decreasing function of the radius R. 

in) E(Vm) is positive for R > Rq close to Ro and converges to —00 when R goes to +00. 

As a consequence, there exists a radius R\ > Rq such that E(Vm) is positive for Ro < 

R < R\ and negative for R> R±. 
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The proof of Prop osit ion 14 . 1 1 relies on technics which are more general than the frame- 
work of this paper. However, to obtain Proposition 14.11 we have to use properties of the 
phase plane corresponding to (12. 3p . These particular properties are observed on the phase 
plane and on the numerical simulations. 

4.1 Energy of the unstable profiles 

We prove here the first assertion of Proposition 14. 1[ The proof relies on the following 
lemma. 

Lemma 4.2. Let V be an unstable profile. Let // > be the first eigenvalue of Ly and 
let <p be a positive associated eigenf unction, the existence of which is stated in Proposition 
\3.b\ Then, there exist two globally bounded solutions U-(t) and u + {t) of ( 12. 4p with the 
following prescribed asymptotic behaviour in ^{—R.R): 

u±{t) = V ±e flt (p + o(e> J ' t ) whent — ► -oo . (4.1) 

As a consequence, for all (y,t) G (—R,R) x M., u_(y,t) < V(y) < u + (y,t). 

Proof: The lemma is a classical application of the theory of stable and unstable manifolds 
near an equilibrium point. We refer to [5] and [2]. Indeed, we can split the spectrum of the 
linearisation Ly into {/i} and the spectrum contained in the half-plane {z G C, Ke(z) < 
\x — e} with e > small enough. As shown in Proposition 13. 6[ \x is simple and admits a 
positive eigenfunction ip. We defined the strongly unstable set by 

W UU (V) = {u G H\-R,R), 3 u(t) global solution of (J22D such that 
lim e(-" +£ / 2 )*|K*)-T/||tfi ( _^)=0}. 

The theory of invariant manifolds near an equilibrium point shows that W UU (V) is an 
invariant one-dimensional manifold which is tangent at V to the line V + Mip. The 
manifold W UU (V) consists in V and two globally defined trajectories u + and w_ satisfying 
the asymptotic behaviour (14.11) . Of course, the last assertion is a direct consequence of 
this asymptotic behaviour, of the positivity of ip and of Lemma 12. 1[ □ 



The first assertion of Proposition 14.11 is deduced from Lemma 14.21 as follows. Let V 
be an unstable profile. We know from Section [3] that V lies between and the largest 
profile Vm, both being stable. Let w_(t) be the solution given by Lemma I4T21 for the 
profile V. By Lemma |2T2"1 the u— limit set of w_ is non-empty and consists in profiles. We 
wonder which profiles may belong to this u— limit set. Assume that there is an unstable 
profile V in u(uJ). This profile must satisfy V < V since U-(t) < V for all t. Moreover, 
by Lemma I2.3[ V ^ V and thus, by uniqueness of ODE solutions, V(y) < V(y) for all 
!/Gl. Let u+(t) be the solution given by Lemma [4.21 for the profile V. For t close to 
— oo, u+(t) < w_(t). Using Lemma [2TT| V < u+(t) < u_(t) < V for all times. Since 
Lemma [2.31 prevents «+(£) to go back to V, U-(t) cannot converge to V, which yields a 
contradiction. Therefore, since w(w_) is non-empty, the only possibility is w(w_) = {0}. 
Lemma O then shows that £{V) > £ (0) = 0. 
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By using u+(t), we prove similarly that £(V) > £(Vm)- Notice that the above argu- 
ments also prove that every unstable profile intersects the other unstable ones as illustrated 
in Figure HI 

4.2 Energy of the stable profile 

The proof of the second and third assertions of Proposition 14. 1 1 is splitted in several parts. 
We use in this section the notations of Section 13.11 with obvious changes to add the de- 
pendence of the different objects with respect to R. 



• i?4 Va/(0) is increasing. 

To study the function R H- V^, we invoke the graphic interpretation of the profiles 
developped in Section 13.11 By symmetry of the largest profile, V^f(O) is the intersec- 
tion of the trajectory y \-t ^ y (V^(—R),cnV^(—R)) with the horizontal axis. There- 
fore, R H- V^(0) is increasing if and only if R \- > Vj^(-R) is increasing. For R 
given, &r(Vm(—R),0!Vm(—R)) belongs to V and at this point of the phase plane, the 
flow $ y crosses T>' transversally from the top to the bottom. Thus, for small e > 
0, $# +e (V^(— R), cxVm(— R)) is below V and by continuity, there exists an intersec- 
tion of § R+£ T> and V at the right of ^ R+£ (V^(-R) : aV^(-R)). This means that 

yR+e 



R — e) > Vm(—R) and thus R >->■ V^f(0) is increasing. See Figure CD 




Figure 7: when R = R Q , a saddle-node bifurcation creates a stable profile V^ (the largest 
profile) and an unstable profile V R . When R > Rq increases, the largest profile V^ is 
increasing. When R — > +oo the largest profile V$ piles up to the singularity (1,0) of the 
phase plane. 

By a standard implicit functions argument, one shows that R H- Vj$ is continuous and 
differentiable as soon as QrD intersects T>' transversally at (V^(R),aV^(R)). We admit 
that it is always the case when R > Rq. 
• E(Vm) converges to — oo when R goes to +oo. 
This fact has already been proved in [3] by the first author. We only give here brief 



11 



arguments. When R increases, the trajectory y G (—R,R) (->• $ y {y R I (—R) ) aV R I (—R)) 
goes to the right and V^(0) converges to 1. For large R, the trajectory y G (—R,R) >->• 
$j/(V^(— R), aV^(— R)) spends a time O(R) in a neighborhood of (1,0), where (u,v) i-> 
\\v\ 2 — F(u) is negative, and a bounded time outside. Therefore E{V^) = S{V^) goes to 
— oo when R goes to +00. 

• There exists (3 G (a, 1) such that R \-> E R (V R j) is increasing if V^(0) < /3 and 
decreasing if V^(0) > (3 

We admit that for all R > Rq, R H- V^(y) is differentiable. Let E p be the energy 
functional defined by (12. 2 p with R = p. Let V be a given fixed profile. We consider the 
variation of p f-> i£p(V) sX p = R. 

d P E p {V)\ p=R = (-F(V(R)) - I |V(i?)| 2 ) + (-F(V(-J2)) - | |f(-.R)| : 

Since F (±i2) = ±y/a V'(±R) and since V"'(y) + f(V(y)) = 0, y ^ \ \d y V{y) | 2 + F(V(y)) 
is constant and 

0,J5 P (V) = -2F(V(0)) . 

Since 1/ is an equilibrium of (12. ip . it is a critical point for the energy £"#. In other words, 
for any W G H^R), 

d h E R (V + hW)=0. 

Henceforth, if R 1— > V$ is differentiable, then combining both derivatives yields 

d R {E R (V*)) = -2F(V*(0)) • 

It remains to notice that R i-> V^(0) is increasing and converges to 1 when i? goes to 
+00. Moreover, there exists /3 > such that F is negative on (0, /3) and positive on (J3, 1]. 
In fact, (/3, 0) is the point of the phase plane where the homoclinic orbit to intersects 
the horizontal axis, see Figure [71 

• E R (V^j) is positive for R close to R . 

When R passes the value Rq, the creation of both non trivial profiles occurs through 
a saddle-node bifurcation: two equilibrium states appear at the same point, one stable 
and one unstable. The stable profile is V^, the largest one, and we denote by V R the 
unstable profile which lies between and V$- From now on, we work in (-R, R) by 
using the analogy presented in Section [2j We want to show that S{V^) is positive for R 
close to R . When R decreases to R : the profiles V$ and V R collide. As £(V R ) > 
for R > R as shown in Section [4.14 we must have at the limit £(V R °) = £(V M °) > 0. 
Assume that £(V^°) = 0. Let R n > R be a sequence of thickness decreasing to R and 
let V Rn be the associated unstable profiles. Let ut{t) be the sequence of solutions of 
( 12. 4 p given by Lemma 14.21 applied to V Rn . We know by Section 14.11 that w™ (t) converges 
to when t goes to +00. Let K = lim n ^ +00 H^^Hh 1 (-#,#)• Notice that K is positive 
since for each n there exists y G (—R,R) such that V Rn {y) > a. We set t n to be a 
time such that \\u"i(t n )\\ H i^_ R ^ = K/2. By a compacity argument similar to Lemma 
12.21 (R is moving, but nothing singular happens), one can extract a subsequence such 
that u_ (t(p(n)) converges in H l {—R,R) to a function u°°. Notice that by construction, 
u°° is neither nor V Ro . The gradient structure of (12.41) shows that, for all n and t, 
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< S{u^{t)) < S(V Rn ) and thus £(u°°) = 0. Let u°°(t) be the solution of flZ3D for 
R = Rq, with initial data u°°. This solution is not a profile, and so its energy S(u(t)) 
decreases and is negative for t > 0. However, u°°(t) must converge to a profile when 
t goes to +oo due to Lemma 12. 31 Since all the profiles at R = R have an energy 
equal to by assumption, this is impossible and we get a contradiction. Therefore, 
S(V R °) = S(V^°) = E(V^°) must be positive. 

5 Discussion 

In this paper, we have proved that there exist two critical thicknesses R\ > Rq > such 
that if < R < Rq, there is no non-trivial profile solution of (jl.3p . If i?o < R, there exist 
non-trivial profiles. One of them is larger than every other one and is stable, whereas 
every other non-trivial profiles are unstable. Finally, the energy of the unstable profiles 
is always larger than the energy of the stable profiles. If Rq < R < R±, the energy of the 
largest profile is larger than the energy of 0, whereas it is smaller if R > R\. These results 
give us informations on the propagation of travelling fronts solution of equation (11.11) as 
stated in Consequence 11.21 

We recall that in equation ( 11.11) . u represents the depolarization of the brain so if 
u(X) = the brain is normally polarized at the point X, and if u(X) = 1 the brain is 
totally depolarized. A depolarization wave in the brain corresponds to a travelling front 
solution of (II. ip where a non-trivial stable profile invades the zero state. Thus, the above 
mathematical analysis of the stability and energy of the asymptotic profiles of equation 
(II. 3p yields informations on the propagation of depolarization waves in the human brain 
during stroke. If the width of the grey matter is smaller than Ri no depolarization 
wave can propagate through the brain whereas if the width of the grey matter is larger 
than Ri, there exist depolarization waves. This may explain why attempts to observe 
depolarization waves in the human brain have received opposite conclusions in different 
studies [TJ [16j [I7J [23j [21]. In the human brain, the grey matter is particularly thin 
and its thickness may vary a lot. The difficulties to observe depolarization waves in the 
human brain could thus be explained simply by the morphology of the brain. If the initial 
depolarization takes place in a part of the brain where the grey matter is large enough, 
then the depolarization can spread through a part of the grey matter. But if the initial 
depolarization takes place in a part of the brain where the grey matter is very thin then 
no propagation will occur. Hence it is possible to observe depolarization waves in the 
human brain [U [TUl \TT\ , but they will not appear in all the cases [221 121] • Moreover the 
depolarizations would not travel over large distances because they will stop as soon as 
the grey matter becomes too thin. For example, they will stop at the bottom of the large 
sulkus, where the grey matter becomes thiner. This has already be observed in the case 
of the migraine with aura. The aura may be due to a depolarization wave [11] [12] and for 
most of the patients the aura stops at the bottom of the Rolando sulkus. 

This paper also gives informations on how the initial depolarization is erased as il- 
lustrated in Figure [TJ If the grey matter is very thin, then the depolarization is quickly 
absorbed uniformly in the excited area. If the grey matter is a little bit larger, the de- 
polarization occurs in a different way. The depolarized area shrinks progressively while 
the cells stay totally depolarized as long as they are in the middle of the excited area. In 
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this case, the repolarization is due to a travelling front where the normally polarized state 
invades the depolarized state. The neurons in the center of the excited area may stay 
depolarized a very long time, which may cause local damages even if no spreading wave 
is observed. To our knowledge this behaviour has never been observed in experiments. 

In order to verify biologically these results and to validate the model, it would be 
interesting to estimate numerically the values of the thresholds Rq and R\. They depend 
on the values of the parameters of equation (11.11) that are hard to compute due to the 
few possibilities of quantitative mesures in the brain. Moreover these parameter values 
may vary a lot from one person to another and from a species to another. Estimations of 
this parameter values for the rodent and for the human may give us more informations 
on the difficulties to observe depolarization waves in the human brain. Indeed even if the 
mecanisms are the same in the rodent brain and in the human brain, the densities of cells 
are totally different and this must influence the values of the parameters of this model 
and thus of the thresholds. For the moment we are just trying to understand numerically 
the influence of each parameters on these thresholds. 

Acknowledgements: the authors would like to thank Thierry Gallay for several fruitful 
discussions. 
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